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Abstract 

The decays of the rj meson are reviewed in the framework of Chiral Perturbation 
Theory. Particular attention is devoted to the electromagnetic channels involving photons, 
either in the final or intermediate states. 
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Several experimental setups have been recently proposed to study the physics of the 
rj meson. Some are already operative, like the SND detector at the VEPP-2M collider in 
Novosibirsk Q and KLOE at the Daphne collider in Frascati Q. Together with the Wasa 
detector at the CELSIUS storage ring in Uppsala, which should start taking data in early 
2002, they will hopefully collect a copious number of r] mesons. The existence of these 
facilities brings the possibility to study many of its decays. Most of the rj decays detected so 
far 1^, like rj — > 77, rj — > tt~^'it~j 01 rj ^ 7r''77, involve external photon(s) and thus proceed 
unambiguously through electromagnetic (e.m.) interactions. Some channels into leptons, 
r] for instance, proceed through photonic intermediate states and share the e.m. 

origin too. The remaining ones, rj Tr~^TT~TT^ and rj — > 3tt^, are purely hadronic but violate 
G-parity (or isospin symmetry), and have therefore small branching ratios, at the level of the 
e.m. ones. We present here a survey on theoretical predictions for rj decays, focused mainly 
on those involving photons. A complete discussion on r/ ^ Svr can be found in this Handbook 
[Q, |5|, In addition, weak decays, with emphasis on C and/or CP violation, are discussed 
in Ref.§. 

Being the rj mass just above the kaon mass, rriri ~ 550 MeV, the variety of e.m. 
r] decays represents, in principle, an excellent laboratory to test Chiral Perturbation The- 
ory (ChPT) 1^], both in the odd and the even intrinsic parity sectors. However, two main 
drawbacks must be mentioned. One is the rj — rj' mixing which requires the inclusion of the 
non Goldstone boson r]i (a pure and massive SU{3) singlet) in the ChPT framework. This 
has been discussed in Ref.Q from the theoretical point of view and no definite predictions 
can be done since, at one-loop, there are many unknown free parameters. The second diffi- 
culty is the presence of counter-terms at next-to-leading ChPT order. Most of the existing 
theoretical predictions assume resonance saturation of the counter-terms. Finally, for some 
processes, definite ChPT predictions are not available and we present instead estimates based 
on phenomenological Vector Meson models. 

We give in Section 1 a brief review on the treatment of the r/-particle in ChPT. Section 
2 is devoted to the odd intrinsic parity rj decays (i.e. those of the anomalous sector, proceeding 
through the Wess Zumino Lagrangian), whereas the even-parity ones are presented in Section 
3. One multi-body odd-parity final state, rj — > ittt^^, is strongly related to 77 ^ vrTrvr and its 
discussion is relegated to Section 4. 

1 77— particle in Chiral Perturbation Theory 

ChPT is discussed in detail by Gasser [Q] and Bijnens |p in this Handbook. Here we give a 
brief theoretical description in order to fix our notation and for self-consistency. 

ChPT provides an accurate description of the strong and electroweak interactions 
of pseudoscalar mesons P at low energies In this sense and at lowest order, ChPT is 
essentially equivalent to Current Algebra (CA) ||lO|] . At higher orders, loop effects appear 
giving rise to both finite corrections and divergences. The former are known to improve the 
lowest order (or CA) predictions, while the latter require the introduction of sets of counter- 
terms -whose number depends on the perturbative order one is working- and restore the 
renormalizability of the theory, order by order in the perturbative expansion. 

The nonet of pseudoscalar mesons P is described in terms of the SU (3) octet and 
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which appear in the ChPT Lagrangian through the parametrization 
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with / = 130.7 MeV 
corrections included [| 



, as can easily be deduced from the charged pion decay, with e.m. 
|3|. The transformation properties of S under chiral ^7(3)^ x U{3)ji 
are given by S — > Ul^uI^ and the physical rj and rj' particles deviate from being a pure SU{2>)~ 
octet and S'[/(3)-singlet, respectively. The r] — r]' mixing is conventionally formulated in terms 
of 6 and/or = + arctan\/2, which refer to the octet-singlet and non-strange-strange basis, 
respectively. 
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7/1 = COS (f) (uu + dd) I \f2 — sin ss 
r]i = sin (p (uu + dd) / \f2 + cos </> ss 
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This formulation of the r] — rj' states mixing has been mainly used in the analysis of radiative 
transitions [12, 13, For even intrinsic parity vertices with 4 pseudoscalars, our treatment 
of the r] — rj' mixing will follow the standard ChPT treatment, i.e., the t] singlet effects are 
included in the 0{p^) counter-term (see below). 

We assume that nonet symmetry gives a good description of the whole nonet except 
for the singlet mass which gets an extra term, correcting the ?7a(1) problem. Because of its 
heavy mass we neglect the singlet component in loops and use the octet component with the 
mass according to the Gell-Mann-Okubo formula. The mixing angle (^) is not negligible and 
can have rather large effects in r] decays. This treatment also fits the mass matrix reasonably 
well if the loop corrections to the masses are taken into account 

The lowest order Lagrangian of ChPT — order two in particle four-momenta or 
masses, 0{p'^), and ignoring the extra singlet mass term — is 



where the covariant derivative 



D^ = d^ + ieA^[Q,n 



(4) 



(5) 



contains the photon field and the quark charge matrix Q [Q = diag(2/3, — 1/3, — 1/3)] 
thus generating the e.m. couplings of charged P's through the commutator The 



*We use the PDG's normalization 
F = is used. 



for the pion decay constant. Notice that in most theoretical papers, 
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non-derivative terms in Eq. (|^ with x 

[M. = diag(m„, m(i, m^)] and lead to 
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2 rriu + ms rriu + ma niu - ma md - mu 
where ?ti^o^ is the mass matrix element describing tt^-t] mixing and 

Ami = i^lo - m\+)QCD = 6.0 x lO^^ Q^y^ (7) 

is the kaon mass squared difference (e.m. self energies removed). Its numerical value is an 
average between the result |^ Ami = (m^o — ''^1+ ~ "^^o + = 5.3 x 10"'^ GeV^, 

following from Dashen's theorem ||l^, and independent estimates []T5|| (including improved 
versions of Dashen's theorem |l^) leading to Ami ™- range (6.5 — 7.0) x 10~^ GeV^. 

The next order Lagrangian, 0{p'^), contains a series of ten counter-terms identified 
and studied by Gasser and Leutwyler Q, 
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and the Wess-Zumino (WZ) term ||T^, responsible for the anomaly and odd intrinsic parity 
processes, whose action (ignoring external fields other than photons) is ||l^, |l8[ 
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In the above equations, F^j^y is the e.m. tensor and the number of colors. 

Accurate numerical predictions require the knowledge of the counter-term values. One 
can rely on experimental data to fix the finite part of their coefficients, as was originally done 
by Gasser and Leutwyler [0] or, alternatively, saturate them by the contributions of low- 



lying resonances (vector mesons usually providing the largest contributions) at ~ Mp |19|. 
Vector mesons can be introduced using a hidden symmetry formulation and assuming ideal 



mixing and nonet symmetry. The vector meson nonet can be described by = p'^\°'/\/2 + 



p^/\/3. The relevant Lagrangian with vector mesons at 0{p^) is 



(10) 
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where cr is a 3 x 3 unitary matrix parametrized in terms of a set of unphysical scalar (com- 
pensating) fields that, when gauged away, generates a mass term for the vector mesons. The 
Lagrangian in Eq.(p!o|) also contains the following hermitian operators 

P/iu = d^p,y - d„p^ + ig[p^,py\, Qf, = -^[^{dfi + ieQA,,)^^ - ^^(9^ + ieQ^^)^], 

(11) 

where ^ = exp[(z//)P]. In order to identify the contribution of the vector mesons to the 
finite part of the renormalized amplitudes one has to integrate out the vector nonet fields p 
assuming their masses to be much larger than all the momenta. Explicit expressions can be 
found in Ref. |13]. 



2 Odd intrinsic parity t] decays 

In this section we discuss the most important odd intrinsic parity decays (sometimes called 
anomalous decays) of the r] meson of e.m. origin. We start with r] — > 77, which gives 
information on the r] — r]' mixing angle, and the closely related ij — > and r] p^ p~ , 

involving leptons in the final state. We consider also the process r] Tr~^TT~j and postpone 
the discussion of the multi-body channel rj mr^yy to the Section 4. 



2.1 ?7 — > 77 and 77* 

Theoretically, -P77 transitions involving on-mass-shell photons, k'^ = k*"^ = 0, contain valu- 
able information on the mixing (quark-content) of the r/, r/' mesons. The amplitudes for 
^8,1 77 in ChPT are 

A{P ^ 77*) = -^^^" gM^^/^g fc^^fc; (12) 

where C^o = 1, C^g = l/\/3 and C^^ = 2\/2/\/3 at lowest order and e and e* are the 
polarization vectors of the photons with momenta k and k* , respectively. 
The present experimental values are Q 

r(7r° ^ 77) = 7.74 ± 0.56 eV 

^ 77) = 0.46 ± 0.04 keV (13) 
r(??'^77) = 4.29 ± 0.15 /fceF. 

Notice that the error for r/ ^ 77 is the largest one, about 10%, considerably bigger than the 
~ 3% error of 77' ^ 77. In fact, the value quoted in (^) for the ry — > 77 decay width is the 
average between two types of measurements which are not in agreement. Indeed, photon- 
photon production measurements obtain 0.510 it 0.026 keV, whereas Primakoff production 
gets 0.324 lb 0.046 keV. A better measurement of this partial decay width is thus clearly 
needed. This would allow to know the total r] decay width with better precision, since it is 
computed from the partial decay rate T{rj — > 77) divided by the fitted branching ratio for 
that mode. 
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It is easy to check that Eqs.(|l^) lead to an unambiguous prediction for the vr*^ —>■ 77 
width of 7.74 eV, in complete agreement with data. The situation is different for the t] and 
rj' two photon decays. We have to remember that rj and ry' consist of a mixture of the r/g and 
r]i, with a mixing angle 9, as quoted in Eq.@. Pioneering works in ChPT proceed by fixing 
f^g to the value obtained at one loop, and used the r] and r]' decay widths to fit 6 and /i. 



They obtained ~ I-S/tt, fru — I-I/tt and 9 ~ —19.5° [|T3|. Most of the existing results for 
radiative rj decays make use of the above values for the mixing angle and decay constants. 

Notice however that the description of the r] — r]' system in terms of a universal octet- 
singlet angle 9 has been shown to be insufficient when a similar pattern is assumed for the 
decay constants. (See Kroll's contribution to this workshop |2^.) Defining decay constants 
through matrix elements of octet and singlet axial- vector currents, one is lead to a description 
involving two decay constants, /^g and /r,^, and two angles 9s and 9i [|^]. They are related 
through 

({I ii)-i'r'; "/'^'"/')- (") 

V/v fri'J v/ssm 98 /icos 9i J 

Alternatively, one can define decay constants related to axial-vector currents with non-strange 
and strange quark content. In this context, one writes 

f?i' f'r,' ) \fq sin (pq fs COS (ps J ' 

Recent phenomenological analysis, including new experimental data for the 7^7 and r]'j form 



factors 1 23] and data on the J /ip ^ rj^ and J/ip ^ rj^ transitions, have performed an overall 



fit to data obtaining |2l| 6*8 = -22.2°, 6*1 ~ -9.1°, /s = 1.28/^, /i ~ 1.20/^, in reasonable 
agreement with the theoretical ChPT corresponding results |2^ 9% = — 2O.5°,0i ~ — 4°,/8 = 
1.28/7r,/i ~ 1.25/7^. Similarly, using the non-strange strange basis, fq = (1.07 ± 0.02)/^, = 
(1.34ib0.06)/,r, i?^!? — (ps — (39.3± 1)°. Notice that the fact that (pq ~ (ps allows to use a unique 
angle (j) ~ 40° in the non-strange strange basis, which is in agreement with phenomenological 
analysis p3| and also with a recent determination at Daphne |2 



For k*^ / one defines a slope parameter of the -P77* transition form factor by 
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(16) 



related to the pseudoscalar structure or dimensions. The tree-level amplitude (|12D predicts 
hp = 0. Experimentally, there is a single measurement in the time-like region, using rj — > 
7;U+^^ [26] leading to 



^ = (1.9 ± 0.4) GeV- 



(17) 



In addition, there are also several sets of measurements of the process 77* — > P in electron- 
positron collisions for rather large values of —k*"^. The most recent set, when extrapolated to 
small values using a pole approximation, finds ]23] 



1.67 ±0.02 GeV- 



(18) 



The next to leading order correction in ChPT was computed for on-shell~photons by 
Donoghue et al. ]12] and Bijnens et al. J27] and the expected cancellation of the divergences 
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was explicitly found. However, when one allows for (at least one) off-shell final photon, 
k*'^ 7^ 0, the cancellation no longer occurs. The amplitude reads 

*2n. ^ + (19) 




+ ^r2 + l[F(r^m^) + F(F^rr^|,)]} 



where the first factor corresponds to the tree-level contribution in terms of physical fp, 
A = i + 1 + ln47r — 7 (7 is the Euler constant) contains the divergent terms, and 



F{m^,x)=m^ [1- - ] \ ln^^^^-^= -2m2, x = — . (20) 

V 4/ V X -a/x-I-Vx-4 

For smah A;*^, F{m?,k*'^) ~ — and, in the case of real photons, /c^ = A;*^ = 0, all the 
effects of loops and non-anomalous £4 terms reduce to change / into fp for P = rjs,r]i. In 
order to give a finite result for virtual photons one has to include counter-terms of 0{p^) in 
the anomalous sector, which are of the type [] 



The analysis was done by Bijnens et al. [27| assuming saturation by vector mesons in a hidden 



symmetry formulation. Using 6 = —19.5°, the slope parameter reads |2 

_ 2/,,+/,3 1 -1 ^^rnlml 1 

where the term inside the bracket corresponds to the chiral loops and the other to the counter- 
terms. Using /^g = 1.3/,rj/r?i = I-I/tt ^-iid fixing the scale // = Mp, one obtains ~ 
[0.19 + 1.70] ~ 1.89 Ge^-^^ being the major contribution due to the counter-terms. 

The agreement with the experimental measurement Eg. (p!7|) , performed in the time-like region, 
is very good. However, it differs substantially from the more precise measurement Eg . (|l8|) , 
obtained in the space-like region. It is amusing to observe that SU{2>) breaking effects can 
be considered not only for the decay constants, but also for the renormalization scale /i. In 
so doing, the relevant scale parameter is not Mp but the mean vector-meson mass, defined 
through its sguared as = (DM^ + M"^ + 2M|)/12 = 0.69 GeV'^ (2|]. Using this scale in 
(]2^ ) one gets 6^ = 1.69 GeV~^^ in nice agreement with Eg. (18). 



2.2 7] fi^fi , e+e 

The rj decays into lepton pairs also are of e.m. origin. At present, data exist for the muon 

channel, with a branching ratio BR{rj fx^ /j,^) = T{rj fj,~^ fj,~)/T(r] all) = (5.8 it 0.8) x 
10~^ Normalizing this result to 77 ^ 77, one gets 

r(r/ 77) 



t A complete analysis of the 0{p^) counter-terms in the anomalous sector can be found in Refs.M, poL ]l 
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where the branching ratio [^] BR{ri 77) = 0.393 it 0.0025 has been used. 

The "reduced" ratio (|2^ ) can be expressed in terms of a dimensionless "reduced" 
complex amphtude R{r] fi^fi^), normahzed to the intermediate r] ^ 'j^ amphtude, leading 
to 

Birj - /.V") = 2(3 - f^^f^-)\', (24) 



where (3 = y 1 — Amj^/rn^. The on-shell 77 intermediate state generates the model indepen- 
dent imaginary part of R 

Imi?(7?^^V) = ^ln[^ = -5.47. (25) 



The unitary bound on i?, B > iJ""**^ is then obtained by setting Rei? = in Eg. ([2^). It 
takes the value 

B"™*(?? ^ fi+fi-) = 1.11 X 10^^ (26) 

In terms of 5"™*, the PDG's result |] reads B{r] /i+/i-)/B"™* = 1 + {ReR/ImRf = 
1.3 lb 0.2 and can be used to extract Rei? from experiment 

Rei2(r,^^V") =±(3.0j:?|). (27) 

From the theoretical point of view and assuming the obvious dominance of the two 
photon contribution, the reduced amplitude R{q^) = R{P l^l^) can be written as |31| 



where = rn^p^p^ = mf and k* = q — k, and F is a generic and model-dependent form 
factor, with F{0, 0) = 1 for on-shell photons. The imaginary part of R is finite and model- 
independent. By contrast, its real part contains an a priori divergent 77 loop (if a constant 
F{k'^,k*'^) = 1 form factor is assumed). The cure to this problem is model dependent, 
and proceeds either through the inclusion of non-trivial form factors -which depend on the 
hadronic physics governing the P — > 7*7* transition- or, in a more modern ChPT language, 
the inclusion of local counter-terms to render the result finite. 

Predictions for several form- factor models are known |3^. Let us mention here the 
simple Vector Meson Dominance (VMD) approach, which essentially implies neglecting direct 
P77 and VFj vertices, thus assuming the full dominance of the chain P — > VV — > 77. The 
corresponding form factor is 

^ = ^^^ = M2 (29) 



which, when plugged in (|28|) , leads to [^, |31f| , 

Rei?(77^/i+/x") = -l-3lo:5, (30) 
in the 5C/(3)-symmetric limit My = Mp^^^. This corresponds to 

S(r? ^ ^V") = (LlSt^e) X 10"', (31) 
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in reasonable agreement with the experimental value (^). For completeness, we also quote 
the corresponding results for the r] e~^e~ decay amplitude: ReR{ri e~^e~) = 31.3 ± 
2.0, Imi?(ry ^ e+e") = -21.9 and ^(t? ^ e+e") = (3.04±0.26)5"™* = (1.37±0.12) x 
10-8. 

From the ChPT point of view, Savage et al. have introduced local Pl^l^ counter- 
terms to render the amplitude finite, 

Cc.t. = ^h^^l^l [xiTriQ^^^d^^ - Q'S^Sts) +X2Tr{Q^^ Qd^^ - Qd^^^Q^)\ , (32) 



where I = e or fi, and Q is the electromagnetic quark charge matrix. Using the MS minimal 



subtraction scheme in dimensional regularization, the "reduced" amplitude reads |31] 



Rei?(,2 = = -^ii^^l±^ - ^ + 31n ^ 

4 2 A 



(33) 



4/3 1 + /? 12/3 /? V/5+1 

where the explicit logarithmic dependence on A reflects the ultraviolate divergence of the 
loop, and cancels with the inclusion of the local counter-terms Xi(^) + X2(A). The present 
experimental r] /U^/U" branching ratio requires a counter-term (for A = Mp = 0.77 GeV) 
given by 

Xi(M,) + X2(M,) = | • (34) 



In turn, the less precise vr'^ — > e^e available experimental data |35| translate correspondingly 
into 

Xi(M,)+X2(Mp) = | ^g^+l . (35) 

Notice that the first values are consistent with the existence of a unique counter-term, assum- 
ing lepton universality, while the second ones are not and, consequently, have to be discarded. 
This strongly relates r] — > fJ,^fJ-~ and vr'^ e^e~ in the sense that an accurate measurement 
for one of the two processes allows a good prediction for the other. 

Adopting the resonance saturation hypothesis, one can go one step further and predict 
the value of the finite part of these counter-terms. This amounts to choose Mp^^ = 0.77 GeV 
as both the subtraction point A and the mass My appearing in our Fyv form factor (^9|). 
The previous result for rj fJ-^fJ-^ ( |3l| ) can be taken also as the prediction from ChPT with 
resonance saturation. In this context, it leads to xi(-^p) + X2{Mp) ~ —14, close to the 



experimental values displayed in the first row of (34,^). An alternative method to estimate 
the counter-terms has been proposed recently using Lowest Meson Dominance in large-A'^c 
QCD [^], with the result xi{Mp) + X2{Mp) ~ —8.8 ± 3.6 , in nice agreement with the first 
value in Eq.(p^). 

We conclude that both the experimental result ([2^ ) and reasonable models indicate 
B{r] i^i^ji^) ~ 10^^, which is well in the range of Wasa possibilities. The corresponding 
predictions for B{7] — > e^e^) ~ 10^^ indicate that this channel might still be detected at 
Wasa. 
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2.3 f] 7r+7r"7 

The lowest order amplitude in ChPT for r]{po) — > 7r+(pi)7r^(p2)7(^) can be easily calculated 
from the WZ action (^). One obtains 

where C'^^ = l/\/3, = \/2/\/3 and e is the photon polarization. The reliability of this 
approach can be tested observing that the related ^tttttt coupling can be expressed in terms 
of the coupling constant F^'^ (C^ = 1 in Eq.(|36|)), which at this lowest order is related to the 
7r°77 coupling, F"" (see Eq.(|l|), C^o = 1), by 

Prom the first equality one obtains F^'^ = 9.7 it 0.1 GeV~^, not far from the experimental 
value = 12.9 it 0.9 ± 0.5 GeV~^. The agreement is improved when the extraction of 

the experimental value is corrected taking into account both higher order chiral contributions 
|38| and e.m. effects. One then obtains F^^''^ corrected = lO-^ ± 1-2 Ge^^ ||]. 

The use of the tree-level amplitude, Eq.(^6|), leads to T{ri 7r"'"7r~7) = 35 eV, 



somewhat lower than the experimental value, Texp{il — > 7r"'"7r 7) = 58 it 6 eV and to a 
photon energy spectrum less softer than the observed one |4C]. Notice that use of fr/g = -as 



predicted in ChPT at this level-, has been made in the above calculation, as well as /^^ = /jr, 
assuming nonet symmetry. Next to leading corrections break this U{3) symmetry results and 



have been estimated to imply fr^^/fn = 1-3, fru/fn = 1-1 IH^ as discussed in subsection 2.1 
where a good description of Tr^,r],ri' 77 was achieved using 9 = —19.5°. The introduction 
of the SU (3) breaking in the decay constants requires one-loop contributions. 

The one-loop diagrams contributing to the processes ij, rj' 7r+7r~7 have been calcu- 



lated using dimensional regularization in Ref.[14|. The results can be conveniently written in 



terms of divergent parts, containing A = l/e + l + ln47r — 7, and the function F in Eq.(pO|). 
The amplitudes for the processes relevant for the r] particle turn out to be 



{1 + T^r^ [A(4m^ - Amj, +pl, + k^)- {Ami + ^P?2) In 4 " 1^ 4 



(38) 



3. iz. ^2 3^^ 3 3 
+4F(m|^,A;2)} 



''Notice that this decay is also analyzed in Holstein's contribution to this Handbook 
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{1 + + 2"^^ + ^P?2) - (4m2 + In ^ (39) 

- (2m|, + IpI^) In ^ + |p?2 + |i"K,P?2) + |^("^x,P?2)]}- 
3 3 

valid for both real {k'^ = 0) and virtual {k'^ > 0) photons. 

The divergences appearing in the loop calculation of the amplitudes are eliminated 
redefining the coefficients of the 0{p^) Lagrangian. The terms contributing to the r]s7r~^7r~'-f 
vertex can be easily read from the general expression obtained in Ref. |0|, The 
relevant terms are 



(40) 



where the first and the second terms take care of the divergences proportional to the pseu- 
doscalar masses and invariant momenta, respectively. The terms of the Lagrangian in Eq. 
( |40| ) are not enough to eliminate the divergences for ryi — > 7r"'"7r~7. This is because the calcula- 
tion in Ref. 11,11] was performed for SU{3)l x SU{3)r broken to SU (3) v. When extending 
the symmetry to include the singlet, new terms such as 

ef^'^'^^A^Mix^^ + Sx^)(5^QSt - Q)^do.^^]tr{^dp^^), (41) 

appear contributing only to this singlet part. 

Numerical predictions cannot be automatically made at this point because the finite 
part of the coefficients in Cq are unknown constants. Here we will assume again that at 
/i ~ Alp resonance saturation works and we fix the value of the relevant coefficients by the 
vector meson contributions a la Bando [20|. Then, the complete expressions for the amplitudes 
of the processes we are studying are 



A{rj, 7r+vr-7) = ^j^^^^e^^^^ A,p,Mo,{l + Cops + ^^12), 



(42) 



im p 

where the last pieces come from resonance saturated counter-terms and C^^op^ and C^g^ips 
are the finite part of the loop corrections quoted in Eqs.(^,39). The k"^ term in these Eqs. 
follows from the (phenomenologically reasonable) assumption that there is no direct (contact) 
coupling between vector-mesons and PPP-states and that the -P77 and PV'j vertices proceed 
exclusively through PVV. 
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The predicted decay width for rj tt^tt^j at 0{p^) in this ChPT approach turns 
out to be r(?7 — > 7r+7r~7) = 47 eV 10, thus improving the tree-level prediction, but still 
1.8c7 too low compared to the experimental value. The correction is now dominated by the 
contribution of the 0{p^) terms of the Lagrangian. Setting their coefficients to be zero, the 
prediction for the decay width would be Tioops{r] tt^tt^j) = 27 ± 3 eV, lower than the 
tree-level prediction and very different from the experimental value. Notice that the above 
results correspond to the rj — ij' mixing angle 9 = —19.5°. If one allows small deviations, 
6 = (—19.5 =F 2.5)°, the decay width becomes T{r] tt^tt^j) = 47 it 5 eV. 

Wasa could provide accurate values of this branching ratio as well as for the corre- 
sponding photonic spectrum (details on the later can be found in Ref.[^]). 



3 Even intrinsic parity rj decays 

In this section we move to the consideration of the even intrinsic parity r] decays (sometimes 
called non-anomalous processes). We start with a brief discussion of the hadronic channel 
rj — > Stt, which is an isospin violating process and contains information on nid — ijiu- Then we 
move to e.m. decays. We consider the decay rj — > 7r'^77, which turns to be very interesting 
from the ChPT point of view. The related decay rj — > , important as a possible 

check of C and CP single photon contributions, is also commented. Finally, the channel 
rj 7r^7r~7r'^7, closely related to ^ Svr via photon radiation, is analyzed. 



3.1 ?7^37r 

The decay r/ ^ Svr was originally discussed in ChPT by Gasser and Leutwyler [^l|. It is 
also extensively discussed in Gasser 's Q, Bijnens [^] and Holstein's contributions to this 
Handbook. Here we include a brief review in order to fix notation that will be used below 
when discussing rj — > 7r+7r~7r''7 and rj tttt^^ decays. The r] ^ Stt transition violates G- 
parity. This makes the amplitude proportional to nid — rriu and therefore vanishes in the 
limit of exact isospin symmetry. Using the ChPT Lagrangian (^) one finds the well known 
tree- level amplitude (CA result) 



A{ri^TT+{p+)TT (p-)7r°(po)) 



B{md 



3^3/2 



1 + 3- 



mt 



mt 



(43) 



where s = (p+ +p_)^, sq = {s + t + ^i)/3 = + m^/3 and B{md — niu)!'! = Amj^ is the 
QCD part of the kaonic squared mass difference (see Eqs.(|6|,0)). Gasser and Leutwyler |41| 
used the Dashen's theorem ITni to fix 



which leads to the lowest order rate 



5.3 X 10' 



GeV^, 



(44) 



^dt{v ^ vr+^-vrO) = 66 eV, 
well below the experimental result Q 

pexp(^ ^ TT+TT-TT^) = 271 ± 26 eV. 



(45) 



(46) 
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This prompted the same authors [^] to introduce the one-loop corrections plus the corre- 
sponding counter-terms. The one-loop expressions for the amplitudes read 



^(r7^7r+(p+)^-(p_)^0(po)) 



B{md - rriu) 
3^3/2 



1 + 3 



•So 



— 



+U+V+W 



A{r] 37r°) = A{r] 7r+7r"7r°) + A{ri 7r~7r°7r+) + A{r] 7r°7r+7r" 



(47) 



where U +V contains the contribution of loop diagrams and 0{p^) counter-terms, while 1 + W 
takes into account the r] — rj' mixing effects (and depends essentially on the counter-term). 
For further reference we quote the corresponding values at the center of the Dalitz plot: 
Uq + Vq = 0.39 - 0.03 + 0.l8i and 1 + - 1-15. The one-loop effects in |l| represent an 
enhancement by a factor 2.4 respect to the lowest order prediction 



160 ± 50 eV, 



(48) 



where the error accounts for higher order uncertainties. This result is in qualitative agreement 



with the analysis done in Ref.|42| where it is shown that the unitarity corrections generated 
by TTTT final state interactions are large. On the other hand, the inclusion of e.m. corrections 



do not modify appreciably the above result [}43 |. 

This result can be reanalyzed having in mind that Dashen's theorem should be cor- 
rected jlC] resulting in Am|^ ~ l.SAmj^ljjrp, which represents enhancements of ~ 70% to the 
previous values, thus going in the good direction to explain the experimental rate. 

On the other hand, some discrepancies with the ChPT predictions seem to appear 
when expanding the decay amplitude around the center of the Dalitz plot. Moreover, when 
comparing the decay rates into neutral and charged pions, the situation is also confusing. 
Indeed, in terms of the branching ratio 



T{r]^ 37r° 



r(r? 



(49) 



one gets nree 
PDG is r, 



' exp 



- 1.51, Vi loop 

1.40 ± 0.03 i 



1.43 l^lj, whereas the experimental value coming from the 
favoring the one-loop result, while a recent measurement at 
the CMD-2 detector reports Q r 



1.52 lb 0.04 lb 0.08, which seems to favor the tree-level 
result, although being compatible with the one-loop prediction. 



3.2 7] ^ 7r°77 

The r] — > 7r''77 decay has a long story, nicely reviewed in Ref. j45U , and furnish a stringent 
test of ChPT. It is also discussed by Bijnens Q. At lowest order in ChPT, one predicts a 
vanishing ry — > 7r''77 decay width contrasting with the measured value |Q |^ 

T{r] 7r°77) = (7.1 ib 1.4) x 10^^ T{rj all) = 0.84 ib 0.18 eV. (50) 

Higher orders in ChPT allow to obtain predictions not far from this value thanks to a rather 
complicated interplay among loop corrections and counterterms. 

recent measurement by Crystal Ball of r{ri 7r°77) = 0.38 ± 0.11 eV was reported by Nefkens 
during this workshop. 
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The amplitude for the decay rj^P) {p)j{<li)'y{Q2) can be written in the following 

form 

M = + af )A + ^(6* + br)B, (51) 

i i 

where A and B are the two kinematically allowed amplitudes 
-4 = (ei • e2){qi ■ ^2) - (ei • q2){^2 ■ Qi) 

B = -(ei . e2)iP ■ qi){P ■ q2) - (ei • P)(62 • P){qi ■ 92) ^^^^ 
+ (ei • g2)(e2 • P)iP ■ qi) + (ei • P)ie2 ■ qi){P ■ ^2). 

Their coefficients a*, bj refer to the tree-level and counter-term contributions, whereas af , bf, 
with P = TT or K, refer to the contributions from pion or kaon loops. The sum extends from 
the lowest order terms, « = 2, to higher ones i = 4,6, ... We also define s = (^1 + 92)^ = 2(?i ■q2, 
t = {P — (?2)^) u = {P — gi)^ and s + t + u = rn^ + m^. With these conventions fixed we now 
proceed to review the different contributions to the amplitude M. 

As previously stated, there is no lowest order 0{p^) contribution from £2 to the 
rj 7r*^77 amplitude, i.e. 02 = ^2 = 0) ^ can be immediately seen from Eqs.(^) and (^) with 
[Q, S] =0 for neutral pseudoscalars. 

At next order, O(p^), the tree-level amplitude coming from Eq.(|8|) is also zero, a\ = 
b\ = 0. This is clearly the case for the WZ-term contributing only to anomalous processes 
and also for the counter-terms in Eq.(^). Indeed, photon couplings are generated either by 
covariant derivatives or by F^i, terms like the two terms proportional to Lg (for one photon) 
and to Lio (for two photons) in (^. But the vanishing of [Q, S] when only neutral vr" and 
rj are involved is enough to guarantee a null contribution in all cases |p7| . This absence of 
counter-terms implies that loop contributions at this order cannot be divergent. They contain 
two vertices of the Lagrangian £2 and a (charged) pion or kaon loop. One obtains 64'^ = 
and the finite expressions 



(53) 




2\/2a / 2 1 2 o 2 \ 2 \ 

"4 = ^^^p i 3s - "ir, - ^rn^ - ^f^K ] H{s, m^) 



with (x = s/m?) 



sH{s,m ) = s dz ay — 5 = — - — — m 



Jo w? — szy 2 2x —Jx + Jx — 4 

(54) 



sH{s,m'i') ~ — — 5- for s « m\. 
24m|^ 



The small correction to proportional to Amj^ has been ignored. 

If one further introduces the numerical values / = 130.7 MeV, Amj^ as in Eq.(0) and 
the r] mass, one obtains the loop contributions up to O(p^) to the r] — > ir^^^ width 

r(4)(r^ ^ vr°77) = 7.18 X 10"^ eV (55) 
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more than two orders of magnitude below the measured width (|50|). The pion loop con 



tribution is small due to the G-parity violation in the vertex r/Tr+vr^vr^, which makes the 
amplitude proportional to — m„ or Am|^ as seen in Eg. ([5^) . The contribution of the 
G-parity conserving kaon loops is also small because the loop-integration introduces the 
function H(s,m\) ~ l/24m|^ with a large denominator in our case with >> s. Further 



contributions are required. This was noted independently by Ecker et al., Ref.|49|. 

Corrections of 0{p^) were partially considered in [p8| , |50|] in terms of resonance satu- 
ration for counter-terms including vectors as well as the scalar ao(980) and tensor 02(1320) 
resonances. Such tree-level contributions were argued to dominate over one-loop corrections 
of the same order using G-parity suppression arguments in the case of pion loops, and the 
smallness of the one-loop functions appearing in kaon loops. The corresponding result was 



r(6)(^ ^ 7r°77) ~ 0.18. eV 



(56) 



At 0{p^) more counter-terms appear and a new type of loop-correction becomes 
potentially important. Since the counter-term contributions are assumed to be resonance 
dominated one can simply take the full VMD amplitude, 



MvMD 



P-q2 



mz 



M2 



t 



+ 



P-qi 



mz 



A 



M2 



t 



+ 



M2 



B 



(57) 



as an "all-order" estimate, i.e., at 0{p^) and higher. This represents a substantial increase 
to the rate, namely, 

TvMD{ri ^ 7r%7) = 0.31 eV, (58) 



in reasonable agreement with older VMD estimates [£T[. Analogous ao(980) and 02(1320) 
contributions are expected to slightly modify this result as discussed above. A new type of 
loop effects looks a priori more interesting. Taking two vertices from the anomalous L\yz 
implies a non-anomalous one-loop correction of 0{p^). Indeed, pion loops are no longer 
suppressed because there is no G-parity violation in the innTj and -qTrn^ vertices and moreover 



the kaon loop suppression due to the function H in Eq.(54) does not occur. However, one 
easily obtains 

r(8)(j^ ^ 7r°77) = 4.8 X 10"^ eV (59) 

to be compared with T^'^^ in Eq.(^). The dominant contribution to r*-^-* comes from the 
absorptive part of the amplitude, which is much larger than that of the pion loop in F^*^) as 
expected. 

A convenient way of presenting the ChPT predictions for F(t/ tt^jj) consists in 
summing up all the contributions which are not negligible and present no sign ambiguities, 
i.e., the non-anomalous pion and kaon loops at 0{p^), Eq.(|^), the corresponding loops at 
0{p^) with two WZ vertices, and the "all-order" sum of VMD counter-terms in Eq. (p7|) . This 
implies |48| 



r{r] 7r"77) = 0.42 ± 0.20 eV, 



(60) 



where the error comes essentially from the sign ambiguity in the interference of the ao(980) and 
02(1320) resonances respect to the vector dominated contribution, and the lack of knowledge 
of higher order contributions. This ChPT result is in agreement with related work by Picciotto 
| |52| (performed also in the same ChPT context) and Ng and Peters [^3[ (who estimated the 
dominant VMD contribution). More phenomenological predictions using constituent quark 
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models exist. Good agreement with the experimental result can be reached with constituent 
quark masses of 300 MeV for up and down quarks |5j]. Measures of the photon energy and/or 
the two-photon invariant mass spectra, which have different shapes for different models, would 
help to elucidate the best theoretical description for ry — > n^jj decay. 

3.3 r] 7r°l+r 

The rare decays ry — > tt^ e~ , ir'^ fi^ fi^ have not yet been observed. They can occur at second 
order of the e.m. interaction without C violation, through the chain rj tt^jj ■7T^e~^e~ and 
TT^ — > The imaginary part of the amplitude is given by forcing the two intermediate 

photons to be on-shell and requires good knowledge of the rjn^j'j decay which unfortunately 
is not the case at the moment, as we discussed in the preceding subsection. In addition, the 
real part of the amplitude diverges unless one introduces suitable form factors for the off-shell 
photons that regularizes the loop integration at expenses of making it model dependent. 

Some years ago, Cheng computed the rj Tr^e+e" decay in the limit mg = 
using a VMD model. He found T{ri vr'^e'^e") = 13 fieV. More recently, Ng and Peters 
have analyzed the rj — > TT^fi~^fi~ channel, where the muon mass can not be neglected. Making 
also use of VMD (and also considering the ao(980) exchange in the two photon intermediate 
state) they give approximate analytical expressions for the imaginary part of the r] — > tt^U 
amplitudes and for the energy spectrum of the pion. The following unitarity bounds were 
obtained 

r-^*(^^,rOe+e-),^,_ = l.llO| f^eV 



VMD (gx) 



The inclusion of the ao(980) exchange affects only the muon channel, leading to O.OIqj 
(0.3lo;2) /^^^ for constructive (destructive) interference. They estimated the dispersive part 
of the amplitude by means of an ad hoc cutoff of 1 GeV to be \ReAf^fi/ImAf^fi\ < 0.5, sug- 
gesting that the rate can not be much greater than the unitarity bound. This conclusion 
is consistent with the results obtained by Heiliger and Sehgal using a dispersive treat- 
ment for the similar processes — > 7r^e^e~ and Ki — > 7r^fi^fj.~, iReA^^/ImA^^l < 0.7, 

i?e^^^//m^^^ < 1.25. 
3.4 r] — > 7r+7r^7r'^7 



The r] — > 7r~^7T~TT^j process has been analyzed in the ChPT framework |57, |5^. At lowest 
order, the amplitude is dominated by pure bremsstrahlung and is thus closely related to the 
isospin violating tj tt^tt^tt^ amplitude. At O(p^), apart from the correction to the inner 
bremsstrahlung contribution, a new structure dependent (or direct emission) contribution 
appears, which does not vanish in the isospin limit. This makes the process interesting since 
it carries new information not accessible from rj ^ Sir. Both contributions were estimated in 
[57 1 using an ad hoc parametrization of the tj — > 7r+7r~7r'^ amplitude. It essentially consists 



in approximating the amplitude by its value at the center of the Dalitz plot, using U + V ^ 
Uo + Vo = 0.39 - 0.03 + O.lSi; I + W V^, value that differs from 1 + Wq ~ 1.15 and 
takes into account the phenomenological r] — rj' mixing for radiative processes, and using the 
enhancement of the amplitude through the correction of Dashen's theorem. This analysis 
was performed considering only the isospin conserving kaon loop contribution, which was 
expected to be dominant. Next order contributions were also estimated using VMD saturation 
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hypothesis, in a similar way as in the rj — > 7r*^77 ChPT prediction. The conclusion reached 
was that the inner bremsstrahlung amplitude strongly dominates. Several technical points 
were refined and improved in |58|. Mainly, a better treatment of the Low's theorem for 
the inner bremsstrahlung amplitude and the inclusion of the pion-loops leads to the result 
BR{r] TT+TT-TT^r, E-y > WMeV) = (3.14±0.05) x lQ-^BR{r] vr+vr-vrO). The contribution 
of the structure dependent contribution to the above result turns out to be a tiny ^ 0.02%. 



4 Tj ^ TTTT'j'j decays 



In this section we comment on the anomalous rj — > Tnrj'y decays. They were first discussed 
at tree- level in ChPT in Ref.|59]. The next order contributions to the neutral channel were 



computed in [SO, 61]. Being multi-body decays, several contributions appear. In fact, the 



amplitude for the neutral channel, for instance, can be written as [31| 



A{7] 



7r°7r°77) = A{r] vr°7r°77)^Q_po;e 

+ A{r] 7r\° -fj\^taii + A{ri 7r°7r°77) 
A{r]^ vr°7r°77)y A/D- 



IP/ 



+ 



(62) 
(63) 
(64) 



The first term proceeds through the chain t] — > vr'^vr'^vr'^ with a subsequent anomalous vr'^ 77 
transition. It thus depends on the non-anomalous and isospin breaking rj — > Stt^ amplitude. 



discussed in Section |3^. Actually, the available phase space suggests that a simplified treat- 
ment of the 7] — > Svr'^, along the lines described in the previous Section is possible. The sec- 



ond contribution proceeds through the isospin conserving decay chain t] — > tt^tt^t] — > 7r'^7r''77 
and turns to vanish in the chiral limit. The third term includes one-particle irreducible 
(IPI) one-loop diagrams, with an r] ^ 7r"'"7r~7(7) vertex from the WZ Lagrangian, followed 
by 7r"''7r~(7) tt^t:^^ rescattering. The last VMD contribution represents an estimate of 



higher-order contributions. The analysis of reference [Sl| can be summarized in the diphoton 
mass spectrum of Fig.l. This figure clearly shows that this partial decay width is dominated 
by the pion-pole in rj ^ Stt^. 

In ^] the charged rj —>■ 7r"*'7r~77 channel was also discussed up to 0{p^). The 
interest of this process from the ChPT point of view is masked by the dynamics of the processes 
T] 7r+7r~7 with an additional radiative photon, and r] tt^tt^tt^ with a subsequent vr'^ — > 77 
decay. 
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Figure 1: Diphoton mass spectrum z = rn^^/rn^ for the ry 7r°7r°77 decay. The dashed hue 
is the vr'^-pole contribution. The upper dotted hne corresponds to the 0{p^) result. The lower 
dotted line is the VMD contribution. The solid line corresponds to the total contribution. 
The symbols at the end of the spectrum show the total result using an approximate expression 
valid at values close to z = 1/4. 
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